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Abstract 

Weierstrass's theory is a standard qualitative tool for single de- 
gree of freedom equations, used in classical mechanics and in many 
textbooks. In this note we show how a simple generalization of this 
tool makes it possible to identify some differential equations for which 
compact and even semi-compact travelling solitary waves exist. In 
the framework of continuum mechanics, these differential equations 
correspond to bulk shear waves for a special class of constitutive laws. 

1 Introduction 

A compact wave is a robust solitary wave with a compact support, outside 
of which it vanishes identically. A compacton is to a compact wave what a 
soliton is to a solitary wave with an infinite support; that is, a compacton is 
a compact wave that preserves its shape after colliding with another compact 
wave. Rosenau and Hyman (1993) introduced these concepts over a decade 
ago and a substantial number of differential equations supporting compact 
waves has since been identified and studied. 

From the mathematical point of view, the emergence of such solutions is 
related to the degeneration of the differential equations of motion at certain 
points, and to the corresponding failure of the uniqueness theorem at these. 
Indeed, compact waves are a patchwork made pasting together at degenerate 
points the different possible solutions (unique in between degenerate points), 
hence patching together non-unique solutions; it follows that they are not 
analytic solutions, and in this respect they are substantially different from 
standard soliton solutions. 

Compact waves are weak solutions of differential equations, which are 
continuous - in contrast to shock waves - but have discontinuous derivatives, 
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similarly to acceleration waves. A simple introduction to compact waves and 
to compactons can be found in a recent article by Rosenau (2005). 

The main objection facing compact waves is that the link between an 
adequate partial differential equation and a constitutive law is often tenuous. 
Indeed, the generic affirmation that compact waves emerge from a balance 
between higher nonlinearity and nonlinear dispersion remains vague and es- 
oteric if it is not supported by a clear and rigorous mechanical derivation of 
the right equations. 

Destrade and Saccomandi (2006 a , 2006 b ) recently proposed a general the- 
ory of dispersive nonlinear acoustics and showed how it is possible to derive 
exact equations governing the propagation of compact shear waves in solids 
with an inherent characteristic length. By applying standard asymptotic 
procedures to these exact equations, it is then possible to justify evolutions 
equations which are similar to the compacton factory known as the K(m, n) 
KdV equation. 

From the mechanical point of view, these results were obtained by a 
careful modelling of the dispersive part of the Cauchy stress tensor. The 
corresponding constitutive equations unify and explain in depth several the- 
ories of weakly nonlocal continuum mechanics, such as the a-LANS theory 
of turbulence, or the Rubin, Rosenau, and Gottlieb (2005) theory of inherent 
characteristic length. 

From the mathematical point of view, the emergence of compact waves 
is investigated by a simple modification of Weierstrass's theory (a useful 
qualitative tool for single degree of freedom equations in classical mechanics). 

In this note we show how another class of constitutive assumptions can 
generate the emergence of compact waves and even of semi- compact waves, 
which are travelling solitary waves with a semi- infinite support. 

2 A generalization of Weierstrass's discussion 

One of the most elegant and powerful tools for the qualitative analysis of one- 
dimensional Lagrangian conservative motions is Weierstrass's theory. When 
an integral of energy (in a generalized sense) exists, Weierstrass's theory 
allows us to understand whether the motion of our Lagrangian system is 
periodic or non-periodic, simply by looking at the roots of the potential 
function. 

Let us imagine that our natural Lagrangian system can be described by 
a single holonomic parameter q, say. Furthermore, suppose that it is time- 
independent. The energy integral of such a system is 



T(q,q)-U(q) 



E, 
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where T is the kinetic energy, U is minus the potential energy (note the 
unconventional choice of sign), and the constant E is determined by initial 
conditions. 

Recall that for natural Lagrangian systems, the most general form of the 
kinetic energy is T = a(q)q 2 /2, where a = a(q) is a positive function of q 
only; also, T = only if q = 0. We will thus assume ([T]) is in the form 

\a{q)q 2 - U{q) = E. (2) 

It follows that our motion is confined to those configurations where U(q) + 
E > 0. The special configurations q (say) where U(q) +E = and U'(q) ^ 
are called barriers because they cannot be crossed by the motion of q(t): they 
split the range of possible values for q into allowed and prohibited intervals. 
Points with U(q) + E = and U'(q) = are soft barriers and separate two 
allowed intervals (see below). 

Solving first (j2J) for q, and then differentiating with respect to t, we find 
in turn 

e = 2 u_M + E i , {) _*M [U{) + E] (3) 

a{q) a(q) a 2 (q) 

Hence at a barrier q which is a simple root of U + E = (that is, such that 
U(q) + E = and U'(q) ^ 0), we have q ^ 0. Such a barrier is called an 
inversion point, because the motion reverses its course after reaching it. 
We now separate the variables in (jHJh and integrate to find 




At a (soft) barrier which is a double root of U + E — (that is, such that 
U(q) + E = U'(q) = 0), the integral diverges. Thus a soft barrier is also 
called an asymptotic point because it takes an infinite time to reach it. 



3 Wave propagation 

Now we turn to wave propagation. Consider the case of semi-linear wave 
equations in the unknown field u = u(x,t), 

u tt - c 2 u xx = F(u) , (5) 

where c is a constant and F a nonlinear function of u. Then for travelling 
waves of speed v, i.e. for u in the form 

u(x,t) = <f(z) , z :— x — vt, (6) 



3 



we obtain the second order differential equation 
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c>" - F{}p) = 0. (7) 



Multiplying by ip' and integrating, we find an equation of the same form 
as the first equation in ([3]), where a is the constant v 2 — c 2 and U is the 
anti-derivative of F; note that E is now related to the integration constant. 

Weierstrass's theory tells us that a periodic wave corresponds to the ex- 
istence of two consecutive inversion points; that a pulse solitary wave with 
infinite tails corresponds to the existense of an asymptotic point followed 
by an inversion point; and that a kink solitary wave with infinite tails cor- 
responds to two consecutive asymptotic points; see Peyrard and Dauxois 
(2004) or Kichenassamy and Olver (1992) for similar discussions. 

Consider the case of the following fully non-linear wave equations, 

u tt - c 2 u xx - c 2 Nh {ul) x = F(u) , (8) 
where cnl is a constant. The travelling wave reduction ffih yields 



v 



2 c 2 ) ip" - 4 L Wf}' = F{fp) . (9) 



Here Weierstrass' discussion must be modified, sometimes to dramatic effect. 

Take for instance the degenerate case v = c . In that first integral 

of @ is 

wr - 2 , do) 

a 

where U(<p) = — J F, E is a constant of integration, and a = 3c^ L /2. 

We can still conduct an analysis a la Weierstrass, but we find that the 
fourth-order power above introduces some new features, not present in me- 
chanical conservative Lagrangian systems. Indeed, now the barriers Tp are 
attainable not only when they are simple roots of U + E = 0, but also when 
they are double roots. This is the case because the analogue to (jl]) is here 



and the integral converges for double roots. 

Saccomandi (2004) gives a detailed discussion on this possibility and ex- 
plains its consequences by the failure of the Lipschitz condition, leading to 
a possible lack of uniqueness, and eventually to the possibility of compact 
waves; see Destrade and Saccomandi (2006 a , 2006 6 ) for examples. 

In the present note we study yet another possibility for Weierstrass's 
theory, touched upon by Ferrari and Moscatelli (1997), Rosenau (2000), and 
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Gaeta, Gramchev and Walcher (2006); namely the case of barriers which are 
roots of non-integer order to the equation U + E = 0. 

For a first glimpse at what may happen in this case, we consider in turn 
two examples with roots of fractionary order. 

Example 1 . In the first example, we take E — 0, a — 2, and 

U(q) = VQ (1 " VQ) ■ (12) 
Here we record two barriers: q l = and q 2 = 1. The integral in (jlj) yields 



t = -2 yjy/q-q + arcsin(2 v /g - 1) . (13) 

Clearly, both barriers can be reached in finite times: with our choice of the 
integration constant, these are t = —n/2 for q t , and t = for q 2 . From (jHJ) 
the acceleration is 

At g 2 = 1 the acceleration is not zero, and the motion reverses; at q 1 — 
however, the acceleration blows up! 

Example 2 . In the second example, we take E = 0, a — 2, and 

17(g) = g 4 / 3 (1 - gV3) . ( 15 ) 

Here we record two barriers as well, again q l = and g 2 = 1. We perform 
the integral in (j4j) for t G [— Sv^, 3v2] and solve it explicitly for g as 



«W = ^-^J • (16) 

Clearly again, the barriers are reached in finite times: q x — at t — ±3^/2 
and q 2 — 1 at t — 0. We also find that the acceleration is given by 

<?=^(4-5g3), (17) 

making it clear that the barrier g 2 = 1 is a configuration associated with 
a finite force, whereas the barrier q l = is associated with an equilibrium. 
Furthermore, note that the right hand-side of (1171) does not satisfy the Lip- 
schitz condition at q = q l = 0. This allows non-uniqueness of solution, and 
in fact we can patch together a compact solitary kink: this is equal to one 
for t < 0, then for < t < 3^2 it is accelerated to negative velocity and 
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Figure 1: Example 2. Compact solitary kink wave, patched together using 
fractionary roots in Weierstrass's discussion, with a potential of the form: 
U(q) = ga(l — gl). 

decreases to reach zero at time t = 3^/2, and then stays there for t > 3^/2; 
see Figure [TJ 

These two examples highlight the complexity and the richness of the 
situations arising when considering non-integer barriers. 

Let us now turn to a more general (non- fractionary) case. We take E = 0, 
a = 2, and 

U(q) = q 2r V(q 2 ) , (18) 

where V, and thus U, have a barrier q > {V(q 2 ) = 0). Then by ©2, the 
acceleration is 

q = rq 2r - V(g 2 ) + q 2r+ l V\q 2 ) . (19) 

We should consider several cases, depending on the value of r > 0. 

If 2r < 1, then q blows up at the barrier q — > 0, similarly to the first 
example above. For this case we conclude that sublinear (2r < 1) roots 
correspond to singular points which cannot be reached in any way. 

If 2r = 1 or 2r = 2, we recover the already discussed cases of simple and 
double roots. 

If 2r > 2, the integral in diverges and the barrier at zero cannot be 
reached in a finite time - it is an asymptotic point. 

If 1 < 2r < 2 however, the integral is finite and the barrier zero is reached 
in a finite time. Also, by (fl9l) . q — at that barrier, and we thus arrive at 
an equilibrium configuration with null velocity, a situation which gives rise 
to a predicament: will the motion settle on this equilibrium configuration ad 
infinitum or will it reverse its course? 
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In fact, the Cauchy problem is ill-posed here, and the accompanying 
lack of uniqueness gives us some latitude to patch together compact waves 
(similarly to the second example above) or semi-compact waves, as is seen in 
the next section. 

4 Semi-compact shear strain waves 

Destrade and Saccomandi (2006°, 2006 b ) show that the motion of transverse 
strain waves in nonlinear dispersive solids is governed by the following equa- 
tion, 

(»W) XX + {aW tt ) xx = pW u , (20) 

where x is the direction of propagation, W is the transverse strain, and p 
is the mass density. The constitutive parameters are p = p(W 2 ), the gen- 
eralized shear modulus of nonlinear elasticity, and a, the dispersion param- 
eter of Rosenau et al. (1995). For simplicity here, a is taken constant and 
the strain wave is linearly polarized, travelling with speed v, see Then 
W = W(x-vt) and flU]) leads to 

pW + av 2 W" = pv 2 W. (21) 

Next we write p in the form 

p(W 2 ) =p - (ap /p) [rW^ViW 2 ) + W 2r+1 V'(W 2 )] , (22) 

where p is the ground state shear modulus, r > 1/2, and V is an as yet 
arbitrary function. For bulk waves, v is arbitrary and here we fix it at the 
sonic speed v = a/ Po/p- Then ( 12TI) is exactly of the same form as ( TL9l . 

We are thus entitled to consider the possibility of barriers of non-integer 
order for the wave. For instance, take V in the form 

V (W 2 )= 1 ^(J m -W 2 ) n , (23) 
a 

where 7 > 0, J m > 0, and n > 1 are constants. Integrating (j2"2"|) with this 
choice gives 

^' 2 = 7 £w^(j m -iy 2 )\ (24) 

Then the following change of variable and rescaling of function 

1 

~ 2 z, u>{Z) = W([z{£)]/y/j^, (25) 



Pjr+n-l 

. ft. 
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give the following non-dimensional version of the governing equation, 



2 , ,2r f-i , ,2 



UJ = UJ 



{l-uj 2 ) n . (26) 



Now we briefly discuss whether fl22|) - fl25|) constitutes a reasonable shear 
response for a nonlinear solid. We recall that the shear stress r (say) neces- 
sary to maintain a solid in a static state of finite shear with amount of shear 
K (say) is r = fi(K 2 )K, given here by 



r 

— = K-7 
Ho 



(J m - K 2 ) n - nK^ (J m - K^- 1 . (27) 



First we see that K must not be allowed to go too far beyond ^/Jm for r 
to remain positive. In practice, this means that for a given material we must 
fix J m beyond the maximal shear allowed before its rupture. Note that the 
actual value of J m has no bearing on the existence and characteristics of the 
shear wave because it does not appear in (1261) . 

Second we remark that the graph of t(K) and the graph of [IqK (corre- 
sponding to a material with a linear shear response) cross - independent of 
7 - for K — 0, K — \JrJ„J (r + n), and K = \[7^ n . The slopes of the t(K) 
graph at K = and at K = \fJm are defined (and equal to /io) when 

2r > 1, and n > 2. (28) 

(Of course these slopes are also defined at 2r = 1 and n = 2, but we leave 
those special cases aside because they lead to simple and double roots in 
Weierstrass's discussion, already treated above.) The slope of the t(K) plot 
at K = y/rJ m / (r + n) is always greater than /i . It follows that between 
K = and K = \JrJ m / (r + n), the solid is strain-softening in shear, and that 
between K = \/rJ m /(r + n) and K = \[T^ ll the solid is strain-hardening in 
shear. 

Third we make sure that the shear response is a monotone increasing 
function of K, by taking 7 small enough so that the equation t'(K) = has 
no root in the [0, \/ J m \ interval. Note that the actual value of 7 does not 
affect the non-dimensional equation of motion (|2"B"1) . 

Figure [2] displays some examples of shear stress responses satisfying the 
requirements just evoked. 

Having checked that the shear response is sound, we may now look for 
solutions to the non-dimensional equation (|26j) describing solitary wave so- 
lutions to the original equation f)20p ; these are homoclinic or heteroclinic 
solution to (!2"B"]) . 

Owing to (1281) . the barrier 57 = 1 is necessarily an asymptotic point. 
Taking r > 1 also creates an asymptotic point barrier 57 = 0, leading to 
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Figure 2: Shear stress response obtained by varying the constitutive param- 
eters in (1271) . The plots cross the linear stress shear response (thin straight 



a solitary kink with tails of infinite extend. However, taking 1/2 < r < 1 
gives a barrier reachable in a finite time. The result is a semi- compact wave, 
coming from value 1 at — oo and decreasing to zero, which it reaches in a 
finite time with zero speed and zero acceleration. It may then remain at this 
value zero. For Figure[3]we took the case r = 3/4, n = 3, and chose u(0) = 
to fix the value of the integration constant. 
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